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ON THE DYNAMICS OF A HIGHER-ORDER RATIONAL DIFFERENCE 
EQUATION 


E. M. ELABBASY, M. BARSOUM AND H. S. ALSHAWEE 


ABSTRACT. The main objective of this paper is to study the qualitative behavior for a 
class of nonlinear rational difference equation. We study the local stability, periodicity, 
Oscillation, boundedness, and the global stability for the positive solutions of equation. 
Examples illustrate the importance of the results. 


1. INTRODUCTION 


In this paper, we aim to achieve qualitative study was of some of the behavior and solutions 
in the non-linear of differential equations 
arl y 


11 EE E EE e 
= y bz! p+ cx, 


= T 2 -.% 


where the coefficients a,b and c € (0,00) while k and £ are positive integers. The initial 
conditions #_;,%_j41,...,Zo are arbitrary positive real numbers such that j = — max {k, 4}. 
Consider a € [0,00),y € II, oo) Qualitative analysis of difference equation is not only in- 
teresting in its own right, but it can provide insights into their continuous counterpararts, 
namely, differential equations. 

There is a set of nonlinear difference equations, known as the rational difference equations, 
all of which consists of the ratio of two polynomials in the sequence terms in the same from 
.there has been many work about the global asympototic of solutions of rational difference 
equations [3], [7], [10], [11],[12] and [13]. 

There has been much investigation and interest in difference equations by several authors 
such Amleh [1] has studied the global stability , boundedness and the periodic character of 


solutions of the equation 
Tn—1 





n41 =a + 
Tn 


Camouzis, [2] investigated the periodic character, and global stability of all positive solution 
of the recursive sequece, 
Una 
Tn+1 = Sik — 
Tn 


Hamza and Morsy in [4] investigated the global behavior of the 
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Ln-1 
n41 = Q + y 
Tn 





Owaidy at al [5] investigated local stability of positive solutions of the difference equation 


Y 
—_ Tr-1 
Unt =a + FI š 
n 





Saleh, [9] investigated the periodic character, invariant intervals, oscillation and global sta- 
bility of all positive solution of the recursive sequece, 
n-k 


Tny FA a A 
n 





In the following we present some basic definitions and known results which will be useful 
in our study. 


Definition 1. Consider a difference equation in the form 
(1.2) Poa SE (ea, PaO) 


where F is a continuous function, while k,l € (0,00). Any equilibrium point T of this equation 
is a point that satisfies the condition © = F (%,Z,X). That is, the constant sequence {£n} 
with x, = T for all n > —k > is a solution of that equation. 


Definition 2. Let T € (0,00) be an equilibrium point of Eq.(1.2). As well as we have the 


(i) An equilibrium point 7 of Eq. is said to be locally stable if for every € > 0 there exists 
o > 0 such that, if x_;,..., 1, £o E (0, 00) with |v_; — Z|+...4+|z_1 — Z|+|x£o — T| < 
o, then |x, — Z| < € for all n > —j. 

(ii) An equilibrium point % of Eq.(1.2) is said to be locally asymptotically stable if it 
is locally stable and there exists y > 0 such that, x_;,...,v%-1,%0 E€ (0,00) with 
|z; — T| +... + a — T| + eo — T| < y, then lim — 

(iii) An equilibrium point % of Eq.(1.2) is said to be a global attractor if for every 
L_7,---,L-1, £0 E (0,00) we have lim t= TA 

(iv) The Eq. (1.2) is said to be globally asymptotically stable if it is locally stable and a 
global attractor of the equilibrium point Z. 

(v) Equation (1.2) is said to be unstable, if it is locally stable equilibrium point at 7. 


Definition 3. The sequence {xn} is said to be periodic with period p, if tn+p = Ln for 
n = 0, 1,..., 


Definition 4. Eq. (1.2) is said to be permanent and bounded if there exists numbers m and 
M with0 <m < M < œ such that for any initial conditions £—j,...,£—1, £o € (0,00) there 
exists a positive integer N which depends on these initial conditions such that m < £n < M 
foralln>N. 


HIGHER-ORDER RATIONAL DIFFERENCE EQUATION 3 


Definition 5. A sequence {x,}°°_, is said to be nonoscillatory about the point T if there is 
exists N > —k such that either x, > T for alln > N or tn < T for alln > N. Otherwise 
{£n} is called oscillatory about T. 


Definition 6. The linearized equation of Eq. (1.2) about the equilibrium point © is defined 
by the equation. 








(1.3) Yn+1 = PoYn + PiYn—k + P2yn—e 
oF (2,2) oF azz oF ,z,2) 
BH, \e, 0,2) , = T, L, T ’ — 
Po ae. Pı — P2 — 


The characteristic equation associated with Eq. (1.3) is 
(1.4) p (A) =A — pod — pyr" — po = 0 


Theorem 1.1. [6|Assume that F is a C!—function and let Z be an equilibrium point of 
Eq.(1.2). We can say that the following statements are true 


(i) The equilibrium point 7 it’s called locally asymptotically stable, if all roots of Eq.(1.4) 
lie in the open unit disk |A| < 1. 
(ii) The equilibrium point 7 it’s called unstable, if at least one root of Eq.(1.4) has 
absolute value more than one. 
(iii) The equilibrium point 7 it’s called source, if all roots of Eq.(1.4) have absolute value 


more than one. 
Theorem 1.2. [8] Assume that po, pı and po € R. Then 
(1.5) lpo] + [pil + |[p2| < 1 


is a sufficient condition for the locally stability of Eq.(1.2). 


2. LOCAL STABLE OF THE EQUILIBRIUM POINT 


The equilibrium point % of Eq.(1.1) is the positive solution of the equation, 





T=Q+ = 
= bx? + cz” 
which gives 
a 
2.1 T= 
(2.1) T=a+ — 


Now let f : (0,00)? + (0,00) be a function defined by 


av 


f (u,v, w) BETTA 
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Then, we have 








Of — > —qyacvlun! 
—* Du T Burrow” 
Of _ yav! 
is) ðv (bw? + cu)?’ 
and 
= Yuy! 
(2.4) of £ abpvľw 


Ow (bw? + cu)? 
Theorem 2.1. If 


2ay < ab+ac+a 


then the equilibrium point T = a + z% of eq (1.1) is local stable. 


Proof. From (2.2)-(2.4), we get 








oF TTT) = — u 
aur") E (abtacta)(b+oq © 
| Si ae ay B 
ao apes 22 
and 
Ol eee ax —aby 
a = = P3. 
rae SY (TO 


Thus, the linearized equation associated with Eq. (1.2) about 7, is 


Yn+1 = Poyn + PrYn—k F P2Yn-t- 
It follows by Theorem 1.2 that Eq.(1.1) is localy stable if 


aby 


—acy eile 
| (ab+ac+a)(b+c) 


(ab+ac+a) (b+c) 


ay 
(ab+ac+a) 


<T 














SO, 





acy + aby + acy + aby < (ab+ ac+a) (b+c), 
which is true if 


2ay < ab + ac +a. 


The proof is completed. 
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3. PERIODIC SOLUTIONS OF EQ. (1.1) 


In this part of the research we are studying the possibility of the existence of periodic 
solutions to the eq. (1.1). 


Theorem 3.1. In the all following cases, Equation (1.1) has no positive prime period-two 
solutions: 


(1) If k and £ are both even positive number. 
(2) If k is odd and £ is even positive number. 
(3)If k is even and £ is odd positive number. 


Proof. Case(1) Suppose that there exists a prime period-two solution 














P,Q, P, q, P; Q; 
If k, £ even then £n = n-k = Tne = q 
(3.1) — 
=Q , 
i b+c 
also, 
a 
3.2 = l 
(3.2) —— b+c 
By (3.1) and (3.2), we have 
aU = p=4 


Similarly, we can prove other cases which is omitted here for convenience. Hence, the proof 














is completed. 
The following theorem states the sufficient conditions that the Eq (1.1) has periodic solu- 

tions of prime period two. 

Theorem 3.2. Assume that k and £ are both odd positive integers and y = 1. If 

(3.3) alc — b) > a?b+ ac + 2abe, 

then Eq. (1.1) has prime period two solution. 


Proof. Suppose that there exists a prime period-two solution 


+P, G,P; 9; Ps G+ 
of (1.1). We will prove that condition (3.3) holds. 
We see from (1.1) that if k and @ odd, then %p_¢ = n-k 
y 
a+ a ; 
bp? + cq? 





p= 
and 

2p ets OU 
bq? + cp?’ 
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we have 

(3.4) bp? + cpq = abp + acq + ap, 
and 

(3.5) bq? + cqp = abq + acp + aq. 


By subtracting (3.4) and (3.5), we have 


b (pP — q?) =ab(p-q)+ac(q-p)+a(p-q), 


then, 

(3.6) (+4) = — 

By Combining (3.4) and (3.5), we have 

(3.7) b (p? +°) + 2cpq = (ab + ac + a) (p +q), 
then, 

(3.8) pP +g = (p+ qY — 2pq. 


Form (3.6), (3.7) and (3.8), we get 











p [teete] + 2pq(c— b) = (ab+ac+a) ad 
ac(ab—ac+a) 
— b(c — b) 
we have, 
w+(p+q)utpqa=0 and (p+q)°—4pq > 0, 
then, 





(Baste) dac (ab — ac + a) 
pees e > 0, 


b b(c — b) 
which is true if 
alc — b) > a®b + a*c + 2abc 


Hence, the proof is completed. 
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4. GLOBAL STABILITY 


Theorem 4.1. Then the equilibrium point T of Eq. 1.1 is said to be global stability. 


Proof. We have the next function 





av 
FQ) = Ot T car 
f non-decreasing for v and non-increasing for u, w 
let m = f (M,m, M) and M = f (m, M,m) 
av 
f (u,v, w) =a + EE 
am” 
(S1) m= OS EMIF MT 
aM” 
4.2 M = — 
oe) bm + em?’ 
from (4.1) 
(4.3) M” (b + c) [m — a] = am’, 
from (4.2) 
(4.4) m’ (b+ c) [M — a] = aM”. 


Subtracting Equation (4.3) of (4.4) produces 
(b+ c) [M7 (m — a) — m’ (M — a)| — a(m? — M”) = 0, 


then 
M=m. 











Hence, the proof is completed. 





5. BOUNDEDNESS OF THE SOLUTIONS 


Theorem 5.1. Let {£n} _ max{k,t} be asolution of Eq (1.1), then the following statements 


are true :- 


(1) Assume that a < b and let for some N > 0,¢ny_j41,-..,Un-1,%n € [2,1] are valid, 


then we have F 
(b+c)a — ab 
by ~~ (b+ c)a7 
(2) Assume that a > b and for some N > 0,@n_j41,--.,%n € [1, z] are valid , Then we 


have a 
y y 
ab — (b+c)a 


bto) "= by 
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Proof. (1) If a < b then xy—041,..-,Un-1, tN € [2,1] 





j aX} y 
Tn =Q ’ 
= ba _; + cx, 
then, 
— IT? 
ab? 
£ ———— ar A 
~ ar(b+c) 
and 
P az} y 
Tn+1 =Q ; 
= bx” + cx, 
then, 
ayi 
— (5) 
bte’ 
(b + c) aT! 
Then zi 
7 a 
(b+c)a De ab 
by T "~~ at(b+c) 


Similarly, we can prove other cases which is omitted here for convenience. Hence, the proof 











is completed. 





6. OSCILLATORY SOLUTION 
Theorem 6.1. Eq.(1.1) has an oscilatory solution If k is odd and £ is even and let k < £, 
Proof. First assume that, 


Co Dip ig Pee T and T ee pag lay < T 











so n 
z arl, 
Tı Q 
ba”, + cxl’ 
é 0 
then B 
a ax? 
Ly QT 
bx? + cx? 
and 
> 2 T 
tı>a+ = 7. 
b+c 
So, we have 
ax 
—k+1 
T = H , 





SO, 


HIGHER-ORDER RATIONAL DIFFERENCE EQUATION 


then, 





Secandiy assume that, 


Pp ag) Cig i <T and T-k+1, es LO > T, 








£ ax! y 
Uy Q 
bx”, + erl’ 
t 0 
then, 
<at ae 
Uy Q — an 
bx? + cT 
and 
a = 
zı <a + —— =T 
b+c 
So, we have 
3 
—— aT k+] 
LI=Q T 





SO, 
>at — 
v2 Q 
bx? + cx? 
then, 
>at a T 
t >a — 
b+e 


One camproceed in prove manwer to show that x3 < T and x4 > ¥ and soon. 


proof is completed. 


7. NUMERICAL EXAMPLES 


Hence, the 














In order to clarify the results obtained, we are offering some numerical examples, as follows 


10 ELABBASY, BARSOUM AND ALSHAWEE 








Example 7.1. Fig. 1, shows that Eq. (1.1) has Local stable solutions if a = b = a = £ 
2,c = k = 7 = 1, £o = 12, £1 = 5, £2 = 2,7 = 22 








Fig.1. 


Example 7.2. Fig. 2, shows that Eq. (1.1) has prime period two solutions if l = k = 1, 
a = (1/16), a = c = 2, b = 1, z2 = 3.3, ae = 0.5.(see Table 7.2) 





















































1 | 3.3000 ||| 17| 1.7741 ||| 83| 1.8027||| 49| 1.8032 ||| 65| 1.8032 
2 | 0.5000 ||| 18| 0.1425 ||| 34 | 0.1344 ||| 50 | 0.1343 ||| 66 | 0.1343 
3 | 1.5974 19 1.7857 ||| 35 | 1.8029 ||) 51 | 1.8032 ||| 67| 1.8032 
4 | 0.3832 ||| 20 | 0.1892 ||| 36 | 0.1844 ||) 52 | 0.1843 ||| 68 | 0.1343 
d | 1.4738 ||| 21) 1.7927 || 37| 1.8030 ||| 53 | 1.80382 ||| 69 | 1.8032 
6 | 0.2656 ||) 22| 0.1873 ||| 388 | 0.1844 ||| 54 | 0.1843 ||| 70 | 0.13843 
7 | 1.5826 ||| 23 | 1.7969 ||| 39 | 1.8031 ||| 55 | 1.8032 ||| 71 | 1.8032 
8 | 0.2220 ||| 24 | 0.1361 || 40 | 0.1343 ||| 56 | 0.1343 ||| 72 | 0.1343 
9 | 1.6133] || 25| 1.7994 ||| 41 | 1.8031 ||| 57| 1.8082 ||| 73 | 1.8032 
10 | 0.1912 ||| 26 | 0.1854 ||| 42 | 0.1843 ||| 58 | 0.1843 ||| 74 | 0.1343 
11 | 1.6792 ||| 27| 1.8009 ||| 43 | 1.8032 ||| 59| 1.80382 ||| 75 | 1.8032 
12 | 0.1702 ||| 28 | 0.1849 ||| 44 | 0.1843 ||| 60 | 0.1843 ||| 76) 0.1343 
13 | 1.7253 ||| 29| 1.8018 ||| 45 | 1.8032 ||| 61| 1.8032 ||| 77) 1.8032 
14 | 0.1565 ||| 30 | 0.1847] || 46| 0.1843 ||| 62| 0.1843 ||| 18 | 0.1343 
15 | 1.7553 ||| 31 | 1.8024 ||| 47| 1.8032 ||| 63 | 1.8032 ||| 79 | 1.8032 
16 | 0.1479 ||| 32 | 0.1845 ||| 48 | 0.1843 ||| 64 | 0.1843 ||| 80 | 0.1343 






























































Tabe 7.2 
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Fig.2. 

Example 7.3. Fig. 3, shows that Eq.(1.1) has Oscilatory solution ifa =b = a 

c= k = 1, to = 2, £1 = 2, £ = 2, T = 2.6666666667.(see Table 7.3) 
1 | 2.0000 17| 2.4925 ||| 33 | 2.4858 ||| 49 | 2.4788 ||| 65 | 2.4698 
2 | 2.0000 18 | 2.9032 ||| 34 | 2.9157 ||| 50 | 2.9302 66 | 2.9478 
3 | 2.0000 19 | 2.4916 ||| 35 | 2.4849 ||| 51 | 2.4773 ||| 67| 2.4686 
4 | 2.6667 ||| 20 | 2.9047||] 36 | 2.9174 52 | 2.9322 68 | 2.9497 
5 | 2.5294 21 | 2.4909 ||| 37| 2.4840 ||| 53 | 2.4763 ||| 69 | 2.4674 
6 | 2.9878 ||| 22| 2.9062 ||| 38 | 2.9191 54 | 2.9342 70 | 2.9521 
7 | 2.5528 ||| 23| 2.4901 39 | 2.4831 55 | 2.4753 71 | 2.4662 
8 | 2.9245 ||| 24 | 2.9077 ||| 40 | 2.9209 ||| 56 | 2.9963 72 | 2.9546 
9 | 2.4936 ||| 25 | 2.4892 ||| 41 | 2.4822 ||| 57| 2.4742 73 | 2.4650 
10 | 2.8885 | || 26 | 2.9093 || | 42 | 2.9227 ||| 58 | 2.99384 14 | 2.9571 
11 | 2.4887 ||| 27| 2.4884 43 | 2.4813 ||| 59 | 2.4731 75 | 2.4638 
12 | 2.8958 ||| 28 | 2.9108 ||| 44 | 2.9245 ||| 60 | 2.9405 76 | 2.9598 
13 | 2.4940 ||| 29 | 2.4876 ||| 45 | 2.48038 ||| 61 | 2.4720 77 \ 2.4625 
14 | 2.9018 ||| 80 | 2.9124 46 | 2.9268 ||| 62 | 2.9427)||| 78 | 2.9624 
15 | 2.4939 ||| 31 | 2.4867||| 47] 2.4798 ||| 63 | 2.4709 79 | 2.4612 
16 | 2.9022 ||| 32 | 2.9141 48 | 2.9288 ||| 64 | 2.9450 ||| 80 | 2.9652 



































Table 7.3 
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Fig.3. 





Remark 7.1. Special cases of Equation 1.1 discussed in the |1] when a = c = y = k = 1, 
b=0 and in [2] when a = b = y = 1, c = 0 and in [5] when a = b = k = 1, c = £ = 0 and in 
[9] when a = b = y= 1, c= L= 0. 
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